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'^ ■ Abstract 

In this paper, we are interested in the transition between regimes where either visco-thermal or 
CL|. multiple scattering effects dominate for the propagation of acoustic waves through a 2D regular 

c/3 ' square array of rigid cylinders embedded in air. An extension of the numerical method using 

Q ' Schlomilch series is performed in order to account for visco-thermal losses. Comparison with 

fi ■ experimental data and results from classical homogenization theory allows to study the transition 
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between a low frequency limit (where viscous and thermal effects dominate) and a high frequency 



c/3 

^ ■ 
o ' regime (where multiple scattering effects become predominant). For this particular geometry, a 



large frequency domain where visco-thermal and multiple scattering effects coexist is found. 
PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

A precise description of sound propagation and absorption in phononic crystal for a 
large frequency domain is of great interest in order to study transition between low and high 
frequency. At low frequencies, i.e. when the viscous and thermal skin depths are comparable 
to the crystal characteristic dimension, important dissipation effects occur. These effects 



ii,y,y,y,y,u. simple 



are well described by the long wavelength homogenization theories 
fluid dynamics problems are to be solved: a Stokes unsteady fluid flow problem for the 
velocity, and a Fourier unsteady heat conduction problem for the temperature. Studies 
describing the propagation by the homogenization theory can be found in literature [7|, |8| . 
When the frequency increases enough, such separation no longer is possible and the problem 
transforms into a full wave propagation problem with multiple scattering effects becoming 
predominent. These effects may lead to band gaps (frequency regions where the propagation 



is forbidden). Different existing numerical methods al 



owing to determine these band gaps 



can be found in the literature J9|, [lO|, llll, ll2|, [l3|, [iJ]. For the purpose of studying the 
transition between low and high frequency regime we first extend the numerical method using 
Schlomilch series 15|, [l6|, llll to include visco-thermal effects. As a result, a numerical method 
is obtained which takes into account viscous, thermal and multiple scattering effects. Then, a 
comparison with experimental data (high frequencies) on one hand and with homogenization 
theory results (low frequencies) on the other hand is performed. Finally the transition 
between low and high frequencies is discussed. 

II. THEORY 

We consider sound propagation and absorption through a 2D periodic square arrange- 
ment of rigid cylinders embedded in air for a large frequency domain. This domain includes 
important multiple scattering as well as visco-thermal effects. The three geometries repre- 
sented on figure [1] are considered. The first problem is a row containing an infinite number 
of cylinders periodically spaced according to y axis and insonified by (right and left) ingo- 
ing plane waves. The second one is an arrangement of A^ rows of cylinders with the same 
periodicity according to x and y axis and subject to the same ingoing plane waves. In a 
last problem, we consider the case of an infinite periodic medium called "phononic crystal" . 
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FIG. 1: Geometries under consideration, (a) one row containing an infinite number of cylinders, 
(b) n-rows containing an infinite number of cylinders, and (c) infinite phononic crystal. 

For the cases (a) and (b), we are interested in the reflection and transmission properties of 
the samples, and for the inflnite phononic crystal (c), we study the dispersion relation of 
the Floquet-Boch modes. The solution of the problems (b) and (c) is based on extensions 
of the solution of the problem (a). In order to account for visco-thermal and multiple scat- 
tering effects, two descriptions are introduced. Firstly, we consider an equivalent surface 
admittance surrounding each cylinder (method A). Using this concept, previously described 
by Tournat et al. [l8|, weak dissipation localized in the cylinder boundary layers can be 
considered. Secondly, a complete description allowing for losses (localized or not) in all the 
fluid is presented (method B). Analytically, the two problems are formulated in terms of the 
displacement u and the (excess) temperature r flelds. In the flrst case (method A), these 
fields derive from an acoustic scalar potential (pa-. 



u 






:i) 



which satisfies the 2D Helmholtz equation [19| (convention e *'^*) 



{A + k 
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(2) 



where uj is the angular frequency and Cq is the adiabatic sound speed. The visco-thermal 
effects occuring in the boundary layers are considered in the following boundary condition 
set at cylinder surfaces. 



{dn + ika(3) 0a = 



(3) 



where dn is the normal derivative and (3 is the equivalent surface admittance operator which 
depends notably on fiuid properties (kinematic viscosity and thermal diffusivity). Expression 



of this operator (in the hmit of small boundary layers) can be found in [18|] for the 2D 
cylindrical case. 

In the second case (method B), the displacement and temperature fields can be re- 
lated ([20|) to three potentials : acoustic </>„, entropic (j)h and vorticity ^, 



^ = V (0a + 0?,) + V A ^ 



(4) 



r 



^a 



Cp PoUh 

where /3o is the coefficient of thermal expansion, Tq the temperature, Cp the specific heat at 
constant pressure and u^ the thermal diffusivity. All these parameters depend on the fiuid 



nature. The three potentials (we set ip = (pv^z) obey to the 2D Helmholtz equations ( 20|) 
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(5) 



where z/„ is the kinematic viscosity. The real wavenumber ka is related to acoustic propa- 
gation while the complex wavenumbers k^ and k^ are related to viscous and thermal waves 
respectively. The three potentials are coupled at the cylinders boundaries by the following 
conditions, 

n = 

(6) 
r = 

cancelling displacement and temperature fields. 

Having in view the determination of the Floquet-Bloch modes (in the x direction of the 
phononic cristal), we will in methods A and B expand the potentials in terms of right and 
left going plane waves. 
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where the "pot" index is either a (for the acoustic potential), h (for the entropic poten- 
tial) or V (for the vorticity potential), and the index n determines the possible wavevectors 
through (/Cpof^n = /Cpot) and the condition {kpot,n)yL = 27m which must be satisfied in order 
to respect the y periodicity. The amplitudes of all potentials are merged into coefficients A 



in figure ([T]). To obtain the reflection and transmission properties of one row (case (a)), tlie 
scattering matrix, wliicli relates tlie outgoing waves to tlie ingoing ones, is introduced. 
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[T] [R] 
[R] [T] 






(8) 



wliere [R] and [T] are tlie reflection and transmission matrices. To determine tlieirs elements 
r^°* ^°* and t^°* ^°* , the outgoing potentials are expressed in two different forms, on each side 
of the row for a given ingoing potential (f)ext{r, n' , pot') = e**^p°*' >"'■''. The outgoing potentials 
can be interpreted as the reflected and transmitted flelds by the row, 
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These can also be interpreted in term of the potential flelds scattered by all cylinders be- 
longing to the row. 
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where the coefficients C^°, ^° and -D^°, ^° are independent of the scatterer (their deter- 
mination is performed numerically according to 2l| for instance). At this point, we must 
estimate Schlomilch series of the form ^°l_oo^"^Hm(A;po(|r — ?^|)e*"* ''"''j. For the usual case 
kpot = ka defined in relation (ISl) with ka G M, and due to the slow convergence of Schlomilch 
series, these sums are developed in the form |l5l . Il6l . Il7l | which present a better convergence. 
For the more original cases kpot = kh and kpot = k^ also defined in relation (^ and having 
non zero Im{kh) and Im{kv), the Schlomilch series are naturally convergent. Then an iden- 
tification between the two couples of equations ([9]) and fITOl) lead to the required refiection 
and transmission coefficients. Consequently the refiection and transmission properties of 
one row (case (a)) are entirely determined. 

Now we consider the case (b) of N rows (a) separated by the distance L. The outgoing 



and incoming wave amplitudes are related by a new scattering matrix, 
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where [Rat] and [T^v] are the reflection and transmission matrices for A^ rows. These two 
matrices can be related to [R] and [T] using the following recurrence relations, 



(12) 




[T](^[ld]-[R^_i][R]j [T;v-i] 

[R] + [T]f[y-[RA.-i][R])"'[R;v-i][T] 



where [Id] is the identity matrix. Using these relations, the reflection and transmission 
properties of A^ rows are simply deduced from the two known matrices [R] and [T] . 

The phononic crystal (case (c)) is interpreted as a superposition of an inflnite number of 
rows. Taking an arbitrary row, the scattering matrix concept (relation ([8])) remains valid. 
Moreover, the potentials obey the Bloch condition 22l |. 
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where the ks are the Bloch wavenumbers to determine. The simultaneous use of the two 
relations ([8]) and (fT3l) leads to the following eigenvalue problem. 
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(14) 



where [0] is the zero matrix. As the reflection and transmission matrices [R] and [T] are 
already known, the dispersion relation giving the solution for the Bloch wavenumbers ks 
can be solved numerically. 

III. EXPERIMENTAL AND THEORETICAL RESULTS 

A flrst comparison is performed at high frequencies (acoustic wavelength comparable or 
less than the lattice parameter L, kaL/n > 1) corresponding to a regime governed by im- 
portant multiple scattering effects. The two models described above allow to compute the 
scattering matrix (^ of one row and then, using equation (TT^ . to predict the transmission 
coefficient of a sample containing a flnite number of rows at normal incidence. This trans- 
mission coefficient can be experimentally obtained. Here, two samples of aluminium rods 
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embedded in air in square configuration (as presented in figured] (b)) are used. The radius 
of tlie cylinders is fixed at R=l mm. Tlie lattice parameter of the sample is L=3 mm. The 
two samples contain respectively 3 and 6 rows. Each of them is 100 cylinders wide. The 
experiment is carried out at ultrasonic frequencies using a couple of broadband transducers 
(centered on 100 kHz) with radius 20 mm. The experimental results are reported on figure [21 
When the number of rows increases, the transmission coefficient tends towards zero inside 
frequency domains corresponding to band gaps defined for infinite phononic crystals. We 
have also reported the predictions of the two models A and B, which are undistinguishable, 
and the prediction for the ideal fiuid (i.e. /5 = in ([3])). The band gap locations are well 
described in all cases because they are essentially due to multiple scattering effects only. 
However, outside these bands, an accurate prediction of the transmission coefficients needs 
taking into account the visco-thermal losses. The two visco-thermal descriptions lead to the 
same predictions because the dissipation effects occuring during the propagation are entirely 
localized in the small boundary layers. Discrepancies between the ideal and visco-thermal 
fluid models are solely due to these weak losses. 

A comparison between our computations and the low frequency (i.e. long wavelength) 
model is now performed when the dissipation effects become more important. Previous 
numerical studies based on homogenization theory (limited to leading order terms) and us- 
ing finite element methods to solve the problem of long wavelength sound absorption in 2D 
phononic crystal (figure [T] case (c)) have been performed by Cortis 7|. The principle of these 
long wavelength simulations is to solve the unsteady Stokes equation for velocity and the 
unsteady Fourier equation for temperature. It exploits in consistent manner the fact that 
(to leading order) the velocity field becomes divergence-free at the pore scale (for the pur- 
pose of evaluating viscous effects) and the pressure field becomes uniform (for the purpose 
of evaluating thermal effects). In order to compare these existing homogenization results to 
those given by the two visco-thermal models including multiple scattering, we compare the 
real part of our complex velocity V associated to the phononic crystal with that given by 
Cortis results. This quantity can be expressed from the real part of the Bloch wavenumber 
ks, solution of the eigenvalue problem flT^ . by the relation : V = uj/Re[kB]- The evolution 
of velocity versus reduced frequency {kaL/n) is presented on figure [3] for a first sample of 
0.90 porosity containing 1 mm cylinders embedded in air. Once again, we have reported 
the results for the ideal fiuid (considering the multiple scattering effects only). At the low 
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FIG. 2: Transmission coefficients of two samples containing 3 rows (upper) and 6 rows (lower) 
versus reduce frequency kaL/ir (corresponding to a frequency range going from 20 kHz to 200 kHz). 
Cylinder radius is R=l mm and lattice periodicity is L=3 mm (thus porosity is <I> = 0.65). We 

have reported the theoretical predictions with (— ) and without ( ) dissipation effects and the 

experimental data (+). 



frequency limit, the two visco-thermal models tends towards Cortis results. Consequently 
and for the geometrical dimension considered here, our two models successfully fulfill the 
two asymptotic limits. This is because in the low frequency regime considered here, the dis- 
sipation effects are still localized in small boundary layers. When the frequency increases, 
multiple scattering effects emerge and, as expected, the (leading order terms) classical ho- 
mogenization results are no more able to predict the correct velocity. Here, the real part of 
the velocity decreases and tends towards zero quickly due to the first band gap where the 
propagation is forbidden. At this point, we can precisely study the transition between a long 
wavelentgh visco-thermal description of sound propagation and the emergence of multiple 
scattering effects. The limit of the classical homogenization theory for a simple 2D phononic 
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FIG. 3: Real part of phase velocity versus reduce frequency kaL/ir. Phononic crystal made up 
of 1 mm radius cylinders (porosity is $ = 0.90). We have reported the theoretical predictions 

considering visco-thermal effects (+) obtained by Cortis, multiple scattering ( ), equivalent 

surface admittance concept (A) and visco-thermal as well as multiple scattering effects (— ). The 
two visco-thermal models considering multiple scattering effects can describe the transition between 
low and high frequency regions. Indeed, in this present case, weak dissipation effects are localized 
at cylinder boundary layers. 



crystal of 0.90 porosity case appears near kaL/n = 0.2. The same limit was found for other 
porosities decreasing from 0.90 to 0.50. Description of sound absorption using classical ho- 
mogenization theory is not accurate beyond this frequency limit. In the presented case, the 
dissipation effects are very weak (1% modification on phase velocity) and localized in very 
small boundary layers. Then, the transition description can be performed by the equivalent 
surface admittance concept (method A) as well as the complete description (method B). 

Now, we translate the frequency values and we consider a 1 /im radius cylinders config- 
uration. This scale law allows us to consider the same multiple scattering effects but more 
important dissipation. Associated results are reported on figure H] in the transition domain 
where both visco-thermal and multiple scattering effects are significant. As expected, the 
observed limit for classical homogenization theory, due to multiple scattering emergence, 
remains unchanged. In this second case, our equivalent surface admittance concept cannot 
describe the transition due to important visco-thermal effects (about 20% modification on 
phase velocity). For this geometry, the transition shall be described by considering the 
complete multiple scattering calculations. 

Besides, we know that in the high frequency limit, the velocity V{u!) given by the ho- 
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FIG. 4: Real part of phase velocity versus reduce frequency kaL/ir. Phononic crystal made up of 
1 iim radius cylinders. In this second geometry, the transition can only be described considering 
important visco-thermal effects not localized at cylinder boundary layers. 

mogenization theory reaches an asymptotic value which is related to the tortuosity a^o p| 
(geometrical characteristic of porous media) with the following relation, 

Co 



lim V{uj) 



'ar 



(15) 



This high frequency limit in view of long wavelength theory is just that matched by the 
ideal fluid results, in view of multiple scattering theory, in the low frequency limit (absence 
of multiple scattering effects). Corresponding values are indicated on figures [3] and H] by the 
top horizontal line. The tortuosity values obtained with the ideal fluid model (at the low 

nn 

frequency limit) are in very good agreement with those from literature [7|, 123|. However, 
these values are never reached when considering both visco-thermal and multiple scattering 
effects as illustrated on figures [3] and HI This observation confirms the existence of a large 
frequency domain (from kaL/ii = 0.2 for this example) where viscous, thermal and multiple 
scattering effects coexist. In this region, sound propagation and absorption can only be 
described by taking into account all these effects together. 

IV. CONCLUSION 

In summary, we have extended an existing model describing the multiple scattering effects 
on 2D phononic crystal (using Schlomilch series) by including viscous and thermal effects. 
The main results consists in the existence of a large frequency domain where both multiple 
scattering and visco-thermal effects shall be taken into account. For an usual millimetric 
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geometry surrounded by air, dissipation effects are very weak in the transition domain. Thus 
the description can be performed using the simple equivalent surface admittance concept 
considered in [l8| (including the visco-thermal effects in the small boundary layers near the 
cylinders). However, for a micrometric geometry, this simple description becomes inaccu- 
rate due to the fact that the boundary layers are not small. In this paper, this case was 
successfully described by a general complete multiple scattering calculation. However, the 
studied transition will also be reproduced by an extended addmittance approach (removing 
the assumption of small boundary layers) as long as the viscous and thermal boundary layers 
are not overlapping between the different cylinders. Furthermore, we have shown that we 
can precisely recover the high frequency limit of the classical homogenization results with 
our complete multiple scattering calculation in the given 2D phononic crystal. 
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